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GENERAL CONVOLUTION IDENTITIES FOR BERNOULLI AND 

EULER POLYNOMIALS 

KARL DILCHER AND CHRISTOPHE VIGNAT 


Abstract. Using general identities for difference operators, as well as a tech¬ 
nique of symbolic computation and tools from probability theory, we derive 
very general kth order {k > 2) convolution identities for Bernoulli and Euler 
polynomials. This is achieved by use of an elementary result on uniformly 
distributed random variables. These identities depend on k positive real pa¬ 
rameters, and as special cases we obtain numerous known and new identities 
for these polynomials. In particular we show that the well-known identities 
of Miki and Matiyasevich for Bernoulli numbers are special cases of the same 
general formula. 


1. Introduction 

The Bernoulli and Euler numbers and polynomials have been studied exten¬ 
sively over the last two centuries, both for their numerous important applications 
in number theory, combinatorics, numerical analysis and other areas of pure and 
applied mathematics, and for their rich structures as interesting objects in their own 
right. The Bernoulli numbers i?„, n = 0,1, 2,..., can be defined by the exponential 
generating function 

oo „ 

(1-1) ^ = (kl<2T). 

— 1 n\ 

n—Q 

They are rational numbers, the first few being 1, —i, i, 0, 0, with 

B 2 k+i = 0 for fc > 1. For the most important properties see, for instance, [TJ 
Ch. 23] or its successor [TH Ch. 24]. Other good references are [TT], [Tl], or [T^ . 
For a general bibliography, see [5]. 

Numerous linear and nonlinear recurrence relations for these numbers are known, 
and such relations also exist for the Bernoulli polynomials and for Euler numbers and 
polynomials which will be defined later. This paper deals with nonlinear recurrence 
relations, the prototype of which is Euler’s well-known identity 

(1.2) ^ BjBn-j = -nBn-i - (n - 1)B„ (n > 1). 
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This can also be seen as a convolution identity. Two different types of convolution 
identities were discovered more recently, namely 


(1.3) 


2 ? D D 

V ^ JD j JDn—j 




n—2 

E 

J=2 


BjBn-j 


Br, 


, .. . = (n > 4) 

jJj{n-j) n 


by Miki where Hn = + - + the nth harmonic number, and 

n—2 n—2 

(1.4) (n + 2)^i?,i?„_,-2^ 

i=2 j=2 

by Matiyasevich [16]; see also [2] and the references therein. These two identities, 
which are remarkable in that they combine two different types of convolutions, 
were later extended to Bernoulli polynomials by Gessel [10] and by Pan and Sun 
[20] . respectively. Gessel [10] also extended (11.31) to third-order convolutions, i.e., 
sums of products of three Bernoulli numbers. Later Agoh [5] found different and 
simpler proofs of the polynomial analogues of ([LSD and (El and proved numer¬ 
ous other similar identities involving Bernoulli, Euler, and Genocchi numbers and 
polynomials. Subsequently Agoh and the first author [3] extended the polynomial 
analogue of El to convolution identities of arbitrary order, and did the same for 
Euler polynomials. Meanwhile, following different lines of investigation, Dunne and 
Schubert [Q derived an identity that has both El and El as special cases, and 
Ghu [4] obtained a large number of convolution identities, some of them extending 
(11.31) and (|1.4I) . 

It is the purpose of this paper to contribute to the recent work summarized above 
and to further extend the identities El and (11.41) of Miki and Matiyasevich. In 
Section 2 we state a general result concerning second-order convolutions, and derive 
some consequences. In Section 3 we introduce a symbolic notation with a related 
calculus, and use it to state and prove a very general identity for Bernoulli polyno¬ 
mials. This is then used in Section 4, along with some methods from probability 
theory, to prove a general higher-order convolution identity which gives the main 
result of Section 2 as a special case. In Section 5 we apply most of the methods 
from Sections 3 and 4 to Euler polynomials and again obtain general higher-order 
convolution identities. Finally, in Section 6, we state and prove several further con¬ 
sequences of each of our main theorems. We conclude this paper with some further 
remarks in Section 7. 


7Z -)- 2 

j 


BjB 


n-j = n{n + l)Bn {n > 4) 


2. Identities for Bernoulli polynomials 
The Bernoulli polynomials can be defined by 

( 2 . 1 ) 

3=0 

or equivalently by the generating function 

( 2 . 2 ) i\z\<2n). 

— 1 n\ 

n—O 

For the first few Bernoulli polynomials, see Table 1 in Section 5. They have the 
special values 

(2.3) B„(0) = B„, B„(i) = (2i-”-l)B„, B„(l) = (-1)”B„, 
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(n = 0,1, 2,...), where the first identity is immediate from comparing (12.211 with 
(HD, and the other two follow from easy manipulations of the generating function 
(1^ . We also require the Pochhammer symbol (or rising factorial) {z)k, defined 
for z e C and integers fc > 0 by 

(2.4) {z)k = ^ p = z{z + 1)... {z + k - 1), 

i (z) 

where the right-hand product is valid for k > 1. 

We are now ready to state our first main result, which will be proved later. 


Theorem 1. For integers n > 1 and real numbers a,b > 0 we have 


(2-5) E 


1=0 


n\ {a)i{b)n-i 


Ij {a + b)r. 


Bi{x)Bn-i{x) = E 


1=0 


a{b)i -I- b{a)i 
I) (a -f 6);+i 

ab 


[a b l)(a -f b) 


BiBn-i(x) 

nBn-i{x). 


The remainder of this section will be devoted to deriving a number of conse¬ 
quences of this general identity. First, it is clear by (USD that we get an analogous 
identity for Bernoulli numbers by simply deleting the variable x. The most immedi¬ 
ate special case is obtained by setting a = b = 1. With (1)„ = n! and (2)„ = (n-fl)!, 
some straightforward manipulations involving the binomial coefficients in (12.511 lead 
to the following identity, which was earlier obtained in [3]. 


Corollary 1. For all n > 1 we have 

(2.6) {n + 2)^^ Bi{x)Bn-i{x) = 2E (?+ 2 )t ^n-i(a;). 

When cc = 0, this identity becomes trivial for odd n since one of Bi^Bn-i will 
be zero except in the cases I = 1 and I = n — 1. For even n > 4, however, we have 
the following identity. 

Corollary 2. For all even n >2 we have 

n n , I 9\ 

(2.7) (n-f 2)E^/^n-z = 2E KE 

1^0 1^0 ^ ^ 


This identity, although different in appearance, is equivalent to dEl- For our 
next corollary we need the shifted harmonic numbers which for real a > 0 and 
integers n > 1 are defined by 


( 2 . 8 ) 


Ha^n • — 


n —1 - 

E^- 


Obviously, FLi^n = Hn- As we shall see, the following result can be considered as 
an infinite class of generalizations of Miki’s identity dOD. 


Corollary 3. For real a > 0 and integers n > 1 we have 


n—1 


( 2 . 9 ) y; 


1=0 


(a);(n - ^ - 1) 

(a)n 


Bi{x)Bn-l{x) = E 


1=1 


n\ a{l — 1)! -|- {a)i 


(a) 


i-i-i 


BiBn-lix) 


4 ” -\- Fid ^Bji(^x). 

a -I- 1 
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Proof. The idea of proof is to divide both sides of (12.511 by b and then take the limit 
as & —>■ 0. On the left-hand side we have for 0 < / < n — 1, 

1 {a)i{b)n-i _ ia)iib+ 1)... {b + n - I - 1) ^ (a);(n - ^ - 1)! 
b (a -I- b)n [a -I- 6)„ (a)„ 

as 6 —>■ 0, and on the right-hand side, for 1 < Z < n, 

1 a{b)i + b{a)i _ a{b + 1)... {b + I - 1) + {a)i ^ a{l - 1)! -|- (o); 
b (ad- b)i+i (a + b)i+i (a)i+i 

as b 0. To take care of the terms that were left out in the limits above, we 
note that the term for Z = 0 in the right-hand sum and the term for Z = n in the 
left-hand sum of (12.5|) combine to give 



a + b 
a + b 


B„(x) 


(a)n 
(a -|- b)n 


Bn{x) 


{a + b)n - {a)n 
{a + b)n 


Bn{x), 


and we have 


(a + b)n - (a)n 

b-i-o b{a + b)n 


dx 


{x)r 


= Ha,r. 


where the second equation follows directly from applying the product rule repeat¬ 
edly to the right-hand side of (12.411 . Putting everything together, we get (j2.9p . □ 


As illustrations of Corollary 3 we state the cases a = 1 and a = 2 separately. 
Corollary 4. For integers n > 1 we have 


(2-10) 


n Bi{x) Bn-i{x) 

2 I n — I 

1^1 






I I 


(2.11) (n + 2)Ea + l)i3/(^)%^=E 


1=0 


1=1 


n + 2 
1 + 2 


B, 


(Z -|-Z-|-1)— 


-|- (n -|- l)(n -|- 2) (^—Bn-i{x) + . 


Proof. After some easy manipulations, (12.911 with a = 1 gives 

(2.12) E Bi{x f^-'}f = E (?) Y^n-i{x) + |i?„-i(:r) -f HMx), 

1=0 ^ 1=1 ^ 


and with a = 2, (j2.9p gives (j2.11ll . We now exploit the symmetry on the left-hand 
side of (j2.12ll and rewrite the sum as 


Bn{x) 1 f 1 

n 9 • ^ 


1 


n — I I 


Bn{x) Bi{x) Bn-i{x) 


+ - Bi{x)Bn-l{x) = 


It V—^ 

+ 4 E 


z 


n — I 


1=1 ^ ' 1=1 

Finally we subtract ^Bn{x) from both sides of (12.1211 and note that HnBn{x) then 
becomes iJ„_ii?„(x). □ 


Using a technique that involves generating functions for Stirling numbers and 
Norlund polynomials, Gessel m obtained (12.1011 as a polynomial analogue of Miki’s 
identity dop. When a: = 0, then we also have symmetry in the sum on the right- 
hand side of (|2.10ll : we can therefore use again the identity l/(n — Z)-|-l/Z = 
n/l{n — Z), upon which we easily recover Miki’s identity. 

Some further consequences of Theorem 1 and Corollary 3 will be derived in the 
final section of this paper. 
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3. Symbolic notation and general identities 

1. The use of symbolic notation in dealing with Bernoulli numbers and poly¬ 
nomials goes back to J. Blissard in the 1860s. Subsequently it was used by many 
other authors, among them E. Lucas in the 1870s and 1880s. Later it was put on a 
firm theoretical foundation as part of “the classical umbral calculus”; see, e.g., [5] 


Here we propose and use a system of symbolic notation that is in some respects 
similar to the classical umbral calculus, but is different in that it is related to prob¬ 
ability theory. Also, this system of notation is more specific to Bernoulli numbers 
and polynomials and (later in this paper) Euler numbers and polynomials. 

The basis for our symbolic notation for Bernoulli numbers and polynomials are 
two symbols, B and 14, which are complementary to each other or, as we shall see, 
annihilate each other. First, we define the Bernoulli symbol B by 


(3.1) = (n = 0,l,...) 

so that, for instance, (EU) can be rewritten as 

(3.2) B„(x) = (t + S)". 
Furthermore, with (HI) we have 


(3.3) 


exp (Bz) = ^ 


n—0 


n! 


- 1 


We also require several independent Bernoulli symbols Bi,..., Bk- Independence 
means that if we have any two Bernoulli symbols, say Bi and B 2 , then 

(3.4) = BkBi. 


Second, the uniform symbol U is defined by 

(3.5) f{x + U)= f f{x + u)du. 

Jo 

Here and elsewhere we assume that / is an analytic function for which the objects 
in question exist. From (13.51) we immediately obtain, in analogy to (13.11) . 

(3.6) = ^ (n = 0,l,...), 

n -I- 1 

and using this, we get 

(3.7) exp(^.) = f:^^"§ = ^. 

n—0 

From (13.31) and (13.71) we now deduce 

00 Yi 

exp (z {B + U)) = ^{B + 77)” ^ = 1, 

n\ 

n—0 

which means that B and 14 annihilate each other, i.e., {B +14)^ = 0 for all n ^ 0, 
in the sense that 


(3.8) f{x + B+L4) = f{x), 
or in other words, we have the equivalence 

(3.9) f(x)=g{x+U) g{x) = f{x + B). 
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Finally, we note that (j3.5l) immediately gives, for any m € K, 

(3.10) uf'{x + uU) = f{x + u) — f{x), 

a difference equation that will be used repeatedly. 

2. It is well known that the Bernoulli polynomials are closely related to the 
calculus of finite differences; see, e.g., the classic books [14] or m- It is therefore 
not surprising that methods from difference calculus turn out to be useful in the 
proofs of our main results. Let be the forward difference operator defined by 

(3.11) Aufix) = f{x+ u) - f{x). 

With two (in general) distinct differences ui , U 2 we compute 

Au^Au^fix) = {f{x + U 2 + Ui) - f{x + Ml)) - {f{x + M2) - f{x)) 

= {f(x + Ml + M2) - /(x)) - (/(x + Ml) - /(x)) - (/(x + M2) - /(x)) 
= ^Ui+U2fi^) ^Uif{x) Au^fix), 
which gives the operator identity 

(3.12) Aui+«2 = + ^U2- 

Similarly, one obtains 

(3.13) A^^^.,22+?i3 - Au^Ay^^Au^ Ay^Ay^ F Ay^Ay^ Ay^ Ay^ 

+ Aui + Ay^ + Ay^. 

To generalize these identities, we use the following notation: For a fixed integer 
k >1 and for any subset J C {1,..., fc}, we denote 

(3.14) Aj:=nA„^., 

j&J 

and we let | J| be the cardinality of J. We can now state and prove the following 
simple but important lemma. 

Lemma 1. For any k > 1 and for real numbers mi, ..., m/c we have 

k 

(3.15) ^liiH—hufc = ^J- 

\J\=j 

The case /c = 1 is trivial, and we immediately see that k = 2 and fc = 3 give the 
identities (13.121) and (I3.13L respectively. 

Proof of Lemma 1. This result can be proved by induction on fc in a straightforward 
way. Alternatively, and more formally, we can use the shift operator 

/(x + m) = e“^/(x), 

with the differential operator d = Then we have A„ = — 1, and 

A„,- 1 = E E n - 1) ’ 


and the result follows. 


□ 
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3. We now apply results from the first two parts of this section to obtain a general 
identity for Bernoulli symbols, and thus for Bernoulli numbers and polynomials. 
In what follows, we assume that for a fixed integer k > 1, ui,...,Uk are real 
parameters. To simplify notations, we write, for a subset J C {1,..., fc}, 

(3.16) uj:=Y\_Uj^ J = {1,..., fc} \ J. 

j&J j&J 

The following is, in fact, a restatement of an intermediate result in [3]. 

Lemma 2. Let ui + • • • + Ufe = 1. Then we have 

(3.17) 

^ (t + uiBi H-+ UkBkT = X! X! ( I T ^ {x + Bo + , 

n! n+1 — 7 )! 

where Bq, ... ,Bk are independent Bernoulli symbols. 

Proof. We apply the operator identity (13.151) to the function 

/(t) := -—{x + Bo + uiBi H-h UkBkT^^ . 

(n + 1)! 

Then the left-hand side of (13.151) gives, with (I3.10L 

(3.18) Au,+...+uJ(x) = Aif(x) = f(x -b 1) - f(x) 

= f'(x+U) = ^ (x + uiBi H-h UkBkT > 

n! 

where in the last step we used (13.8p . i.e., 6o is annihilated by U. Similarly, we have 
for any i = 1,... ,k, again using (I3.10L 

(3.19) A„,/(x) = /(x -b Ui) - f{x) = Uif'{x + UiU) 

= Ui— (x + Bo + (uS){i_,,,_fe}\{j})" , 

having used the fact that the uniform symbol U annihilated the Bernoulli symbol 
Bf, note that the coefficients ui have to match for the annihilation (i.e., identity 
(j3.8[l 'l to apply. Using the definition (13.1411 and successively applying ()3.19|1 . we get 

Aj/(x) = -- ' {x + Bo + (uS)7)"+^“^ . 

[n-j + iy. 

Finally, applying (13.151) to this and to (13.181) . we immediately get (13.1711 . □ 


While the case fc = 1 is trivial, for fc = 2 and fc = 3 we get the following identities. 
Examples. For ui -b M 2 = 1, we have 

(x -b UiBi + U2-B2)" = Ui{x + Bq + U2-B2)" -b U2 (x -b So -b MlSi)" 

-b uiU2n (x -b So)"~^, 

and for ui-\- U2 + U2 = 1, 

(x -b uiBi + U2B2 + uaBs)'^ = [ui {x + Bq + U2B2 + U3B3)"' + o.t.] 


nuiU 2 (x -b So -b usBs)'^ ^ + o.t. -b n (n — 1) U 1 U 2 U 3 (x -b So) 

where “o.t.” in each of the first and second rows stands for the “other terms” 
obtained by cyclically permuting the subscripts {1, 2,3}. 


.n-2 
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Remarks. (1) The left-hand side of (j3.17p . and in fact also the terms on the right- 
hand side, could be written as Bernoulli polynomials of higher order, as defined in 
identity (30) in [5] p. 39]. We will not pursue this further. 

(2) It is clear from the proof of Lemma 2 that, more generally, for any analytic 
function / and ui -I- • • • -I- itfe = 1 we have 

n 

f{x + uiBi H-h UkBk) = X! X! {x + Bo + {uB)j ), 

\J\=3 

and in particular, for fc = 2 and ui -I- U 2 = 1, 

f{x -I- uiBi -I- U2B2) = U2f{x + Bo+ uiBi) + uif(x + Bq + U2B2) 

+ UiU2f{x + Bo). 

The main results of this paper are based on Lemma 2 and an analogue for Euler 
polynomials, and will be obtained by considering the expectation when ui,...,Uk 
are taken to be certain random variables. 


4. Generalization and proof of Theorem 1 


1. In this section we prove a higher-order analogue of Theorem 1, of which the 
latter is an immediate consequence. The proof uses some probabilistic methods 
which will be summarized in a brief subsection. 

Theorem 2. For integers k > 2 and n > 0 and for positive real parameters 

Oi,..., Ofc we have 

(4.1) 


E 

/i_i— 


(ai)q ... {ak)ik 
h, ■ ■ ■, IkJ {a-1 -f ■ • • -b Ufe)) 


■Bi,{x)...Bi^{x) = X! E 


ajn\ 


1=11 J|=i 


(n-b 1 - j)! 


E 


-—j 

=n+l-j 


+ 1 ~J ^ (®ij+i)ii ■ • ■ (Qifc)Zfe- 


^ 0 , ^ 1 , ■ • ■ , ^k—jj (ni -b ■ ' * -b ak) 


n+l — lo 


-Big{x)Bi^ ...Bi 


k-j ■ 


When fc = 2, we immediately get Theorem I. For fc = 3 and oi = 02 = 03 = 1 
we get, after some easy transformations and renaming the summation indices. 


(4 


.2) (n + 3) ^ B,ix)Bj{x)Bi{x) = 3 ^ ^ B,{x)B,Bi 

— n d-L-iJ-l — r) ^ ' 


n -b 3 
i 


Bi{x)Bj -b 


n + 3 
5 


Bn-2{x), 


+3 E 

valid for n > 3; this is Corollary I in [3]. Other special cases with k = 3 will 
be considered later, in Section 6 . For arbitrary k > 2, with oi = • • • = ak = 
I, we recover Theorem I in [3], which for a:: = 0 gives a fcth order analogue of 
Matiyasevich’s identity dm, namely 

k 


Y. .= 


li-\ - hlk—n 


n + k 


1=1 


E 


^o+^iH-j 

^n+l-j 


n + k 

lo 


BigBi^ .. .B, 


^k — j 


2. We now summarize some facts from probability theory that will be used in 
the proofs that follow. For the basics we refer the reader to any introductory text 
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in probability theory, e.g., [7] or [23]. For the interplay between probability theory 
and umbral calculus, see [22] , 

We assume that X is a continuous random variable with probability density 
function fx{x), i.e., 

(4.3) Pr(X <x) = f fxi.y)dy. 

J — OO 


Given a measurable function g : M ^ R such that the image random variable g{X) 
is absolutely integrable, its expectation can be expressed as 


(4.4) 



9{y)fx{y)dy. 


The main tool in this section is the use of random variables with a gamma distri¬ 
bution of “scale parameter” 1. We write such a random variable as X ^ Fa, with 
“shape parameter” a > 0, defined by the density function 


(4.5) 


fx{x-,a) 


{ 1 

0 


for a; > 0, 
otherwise. 


Then from the definition of the gamma function, 

^OO 

r(s) = / x^~^e~^dx, 

Jo 

and with (14.4|) and ()4.5|1 we immediately get, for an integer n > 1, 

(4.6) E(r") = y-^y-^e-ydy = = (a)„. 

An essential property of the gamma distribution is additivity, i.e., if Fq^, ... ,Fa,i, 
are independent gamma distributed random variables, then 


(4.7) 


Fa, + • • • + Fq 


F 


oiH- \-ak 5 


where the symbol ^ indicates that the random variables on both sides have the same 
distribution. The relation (HTl) follows from the fact that the density probability 
function for the sum of two independent random variables is the convolution of the 
individual ones; see, e.g., [23l p. 107]. 

The next important tool is the choice of random coefficients ui,... ,Uk such that 
(mi, ..., Uk) follows a Dirichlet distribution with parameters (ai,..., Ofc). This is 
equivalent to choosing k independent gamma random variables Fa^, each having 
shape parameter a^, and to define 

(4.8) Ui = - — , I <i<k; 

^ ai “T ■ ' ' “T -L ttfc 

note that ui + ■ ■ ■ + Uk = 1- For Dirichlet distributions in general, see, e.g., m 
p. 231]. We now need an important property of gamma random variables, namely 
that Fq -|- Ff, and ra/(ra + F;,) are independent when Fq and F;, are. In fact, this 
characterizes gamma random variables; see [15]. This is easily extended to the 
statement that 

(4.9) Fai-f -hFaj, and -— , l<i<k, 

i ai “T ' ■ ' "T -L afc 

are independent. The importance of this lies in the fact that E(Ay) = E(A)E(F) 
for independent random variables X and Y. 
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Combining all of the above, we first note that for any positive integers h,... ,lk 
we have by (14.81) . 

(4.10) E[(r,, + •. • + ... 4'“)] = E(ri^^... 

= E(rig...E(ri';). 

On the other hand, using the independence of the terms in (14.9L we see that the 
left-hand side of (14.101) is equal to 

(4.11) E(ri^X-+'^JE(u^..4'=), 

having also used (14.71) . Finally, applying (14.61) to the right-hand side of (|4.10l) and 
to (I4.11|) . we get 


(4.12) 



(Qi)ii ■■■{ak)h 

(oi -I- • • • -I- a/c)q_|— 


This identity will be used repeatedly in what follows. 


3. We are now ready to prove Theorem 2; as we shall see, much of the work was 
already done in obtaining the identities (13.171) and (I4.12p . 


Proof of Theorem 2. We choose ui,... ,Uk as in (14.81) . Since ui -I- • • • -I- Ufc = 1, we 
can rewrite the nth power term on the left-hand side of (13.171) as follows, and then 
apply a multinomial expansion, using dUl): 


(4.13) 


{ui{x + Bi) -\ -1- Uk{x + Bk)T 


= E 

ll-\ -71 



...u[^Bi^{x)...Bi^ix). 


Similarly, we use multinomial expansions for the powers on the right of (13.171) , this 
time combining the terms a: -I- i3o for the sake of applying (13.21) : 


(4.14) {x + Bo + {uB)j) 


n-j+l 


E / ') Bioix) {u,^+,Bi^+,y\.. {ui^B,^) 


^0+^iH- 

=71+1-j 


— i 


E 


^o+^iH- 

=71+1-j 


n-+l — j \ I Ik-j n I \n u 

I I , Bi,ix)Bi,...Bi,_ 


where we have also used (13.4p . All that remains to be done now is to compute the 
expectation on both sides of (I3.17L which mainly involves applying (14.121) to the 
right-hand sides of (14.131) and (14.141) . In particular, keeping the first notation in 
(I3.16P in mind, we have 


E(ujM 


h 

ij+1 



(oii)i ... (uq )i (aq+i)q ... 

(oi -)-••• + ak)n+l-lo 
(«! ak)n+l-lo ' 


where we have used the fact that (a)i = a and, in the denominator, that Id-.-.-l- 
-I- ■ • ■ -I- Ik-j = j + {n — j + 1) — lo = n + 1 — Iq. This completes the proof. □ 
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5. Euler numbers and polynomials 


The Euler numbers and polynomials are often considered in parallel with their 
Bernoulli analogues. Indeed, they are similar in various respects, including their 
importance in the classical calculus of finite differences (see, e.g., [14] or [H]). In 
this section we follow the outlines of the previous sections to derive analogous results 
for Euler polynomials and, to a lesser extent, Euler and Genocchi numbers. 


1. The Euler numbers En, n = 0,1, 2,..., can be defined by 


(5.1) 




n—0 




The Euler numbers are all integers with = 0 when n is odd; the first few values 
are listed in Table 1. The Euler polynomials can be dehned by 


(5.2) 


Euix) 


3=0 


or equivalently by the generating function 
2e"^^ 


(5.3) 


= ^En{x)- 
n=0 


(|z| < tt). 


A key consequence of (|5.3I) is the functional equation 

(5.4) Er,{x) + Er,{x + 1) = , n = 0,l,2,..., 


which gives rise to numerous applications. One important difference to the Bernoulli 
case is the fact that i?„(0) is not the nth Euler number. The Genocchi numbers 
G„, are often used instead; they are closely related to the Bernoulli numbers via 

(5.5) := 2(1 - 2’^)B„ (n = 0,1, 2,...). 


These numbers are all integers; the first few values are also listed in Table 1. 


n 

Bn 

En 

Gn 

Bn{x) 

En{x) 

0 

1 

1 

0 

1 

1 

1 

-1/2 

0 

1 

X — 2 

1 

H 

2 

1/6 

-1 

-1 

X^ — X + ^ 

1 

H 

3 

0 

0 

0 

x^ — + \x 

„3 32 1 1 

tb 2 4 

4 

-1/30 

5 

1 

x'^ — 2x^ ~ 'io 

x^ — 2x^ + X 

5 

0 

0 

0 

_ —'T 

*Xj 2 1 ^ ^ tXj 

tb 2 "-b 1 2 2 

6 

1/42 

-61 

-3 

^6 _ 1 ^1 

cC KjtXy \ 2 ^ 2 ^ ^ 42 

T® — 3a:^ + 5a:^ — 3a; 


Table 1: Bn, En,Gn, Bn{x) and En{x) for 0 < fc < 6. 


By elementary manipulations of the relevant generating functions, we get 
(5.6) E„(0) = ;^G„+i, E„(i)=2-"E„ (n = 0,1, 2,...). 

The Euler polynomial analogue of Theorem 1 can now be stated as follows. 
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Theorem 3. For integers n > 1 and real numbers a,b > 0 we have 

n 

(5.7) j;: 


1=0 


n\ {a)i{b)n-i p, / 5 p / 5 4 

' Ei\x)En-i{x) = — —Bn+i{x) 


IJ (a + b)r 


n + 1 


n+1 


n+1 

E 

1=0 


+ (a)z + {h)i 

I ) {a + b)i 


Ei{0)Bn+l-l{x) 


As in the case of Theorem 1, this result follows from a higher-order convolution 
identity that will be proved later. As a special case of (15.71) . for a = & = 1, we get 
the following Euler polynomial analogue of Matiyasevich’s identity: 

" /n\ /n 2\ 

(n + 2)^ r\Ei{x)En-i{x) = + 2)Bn+i{x) - i ) Bi{x)En+i-i{0). 

This identity was earlier obtained as Corollary 2 in [3]. 


2 . As we develop a formalism parallel to that involving the Bernoulli symbol, 
we note that the analogue of Bn is £'„(0). Thus, we define the Euler symbol S by 


(5.8) 8^ = En{0), n = 0,1,2,..., 

and elementary manipulation of the generating function ()5.3I1 gives 

(5.9) En(x) = (x + Sr, n = 0,1,2,.... 

The analogue to the uniform symbol It defined in Section 3 is the uniform discrete 
symbol V with generating function 

(5.10) ^ 
or equivalently defined by 

;(, + V)=/(£)±|(£±i) 

for an analytic function /; this is a discrete analogue of (13.51) . With a change of 
variable we have for any real u, 

(5.11) f{x + uV) = ^f{x) + ^f{x + u), 

which is analogous to (13.101) . and which will be just as useful. Next, by multiplying 
(15.31) . setting x = 0, with (I5.10L we see that in analogy with (13.81) and (13.91) we 
have f -b V = 0 in the sense that 


(5.12) fix + £ + V) = f{x), 
or in other words, 

(5.13) f{x)=g{x + V) O g{x) = f{x + £). 

3. The functional equations (15.41) and (15.111) give rise to the definition of the 
discrete forward difference operator 5^ defined by 

(5.14) 5./(x) = + + 

Thus, in particular, we have 5iEn{x) = x" and by (15.lip . 

(5.15) Sufix) = f{x+ uV). 
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In analogy to (13.121) we now compute 

or r ft . fix) + fix + U2) fix + Ul) + fix + U2+Ul) 
2 duJuJix) = - - -+- - - 

_ fix+) + fjx + m + U2) fix) + fjx + Ui) 

~ 2 ”^2 

fix)+ fix + U2) , 

+ - 2 - 

which gives the operator identity 

(5.16) Sui-\-U2 — 2SuiSu2 ^UI ^U2 1- 

Similarly, one obtains 

(5.17) 5 Ui-{-U2+U3 — ^^Ui'^U2^U2 2Su^Su2 26u2^U2 

“t“ ^U\ “t“ ^U2 d” ^U3 • 

In general, using the notation Sj, with the same meaning as in (13.141) . where again 
J C {l,...,fc}, we have the following result. 


Lemma 3. For even k >2 we have 


(5.18) 

and for odd k > 1, 

(5.19) 


= 1 - E E 

1=1 l- 7 |=i 

k 

Su,+...+u, = E E i-^y~"^j- 

1 = 1 1 Ji=j 


These identities 
duction beginning. 


can be proved by straightforward induction, with (I5.16P as in- 
Using notation from (13.161) . we now obtain the following result. 


Lemma 4. Let ui + ■ • ■ + rife = 1. Then for even k > 2 we have 

k 

( 5 . 20 ) in + 1) (a; + Ui£i + • ■ • + Uk£k)^ = E(“2)'^ E ix + B + iu£)j)'^~^^ , 

|J|=i 

and for odd k > 1, 


(5.21) (a; + ui£i H-h Uk£k)'^ = E(“2)'’“^ E (a; + + iu£)j)'^ , 

a=i \J\=3 

where £q, ... ,£k are independent Euler symbols. 


Proof. Comparing (15.141) with (13.111) . we get the operator identity A„ = 2i5„ — 2, 
and thus for even k > 2 we have with (I5.18p . 

k 

(5.22) A„j-|— \-u^ = E](“2)-^ E] 

1=1 \J\=J 

Since ui + ... Ufc = 1, we have by 1)3.18|) . 

(5.23) A„,+...+„,/(ai) = /'(x+^^). 


Now let 


n+1 


fix) := ix + B + Ui£x H-h Uk£k) 
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and apply (|5.22ll to this function. On the left-hand side, using (j5.23|) . the symbols 
U and B cancel each other, and we get the left-hand side of (j5.20|) . To obtain the 
right-hand side, we first note that for any i = 1,..., fc we have by (j5.15l) . 

(5-24) SuJix) = {x + B+ , 

having used the notation in (I3.16|) and the fact that the symbols V and £i cancel 
each other. As in (13.1911 . the coefficients Ui have to match for this cancellation to 
apply. Successively applying (15.241) and using the notation (I3.14|) . we get 

Sjf{x) = {x + B+ (uf) . 

This, combined with (I5.22L completes the proof of (15.201) . 

The proof of (15.211) is very similar: Instead of (15.2211 use (j5.19l) and apply it to 

f{x) := {x £q + Ui£i -I- • • • -I- Uk£k)'^ ■ 

While the right-hand side is evaluated as before, for the left-hand side we use (15.151) 
with u = 1. □ 


4. We are now ready to state and prove the main result of this section. 


Theorem 4. Let n > 0 and k > be integers and ai,.. .,ak positive parameters. 
Then for even k >2 we have 


(5-25) E (, " 

l,+.^,=n V^l’ ...Jkj (ai + ---+ ak)n 7^1^+^ 

- \-ak)n+i-io 


and for odd k > 1, 


(5-26) E (, " = 


S.i: L..: 




i^ij + i)h ■ ■ ■ i.^ik)lk-i 
' 0 ; • ■ • ) (^1 4 “ * ’ * 4 “ O.k'jn—lQ 




For k = 2, the identity (I5.25P reduces to Theorem 3. In the special case 
oi = • • • = Ofc = 1, the identities (15.251) and (15.261) reduce to Theorems 2 and 3, 
respectively, in [3]. Other special cases can be found in Section 6. 


Proof of Theorem 4- The proof is almost identical to that of Theorem 2: We ex¬ 
pand the powers on both sides of (15.201) and (15.211) using the multinomial theorem, 
and then compute the expectation on both sides by way of (j4.12p . having chosen 
ui,... ,Uk as in (1481) . 

The left-hand sides of (15.2511 and (j5.26ll are obtained just as in the expansion 
(I4.13|l . with Euler instead of Bernoulli symbols and polynomials, and having used 
(15.911 in place of (13.211 . The right-hand sides are expanded as in (14.1411 . with the 
appropriate exponent and with “Bernoulli” replaced by “Euler” where appropriate. 
Applying (14.121) then completes the proofs of both identities. □ 
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6. Some further identities 

In this final section we state and prove some further consequences of our main 
results from Sections 2, 4 and 5, respectively. 

6.1. Consequences of Theorem 1. In addition to the two identities in Corol¬ 
lary 4, we can obtain one more consequence of Corollary 3 by multiplying both 
sides of (12.91) by a and then taking the limit as a —>■ oo. Then all the terms in the 
sum on the left disappear, with the exception of the I = n — \ term. On the right, 
the fraction in the sum tends to 1 for Z > 2, and to 2 for Z = 1. Putting everything 
together, we get the following consequence. 

Corollary 5. For integers n > 1 we have 

(6.1) ^ = n{x - l)B„_i(a;) - (n - l)B„(a;). 

For a; = 0 this is Euler’s identity (II.2L but it is also a special case of identity 
(5.11.2) in [12]. 

We can obtain even more consequences from Theorem 1 by setting a = b = e 
and then taking the limit as £ ^ oo, or by considering the terms in (12.51) as power 
series in e. We begin with the first case. 

Corollary 6. For integers n > 1 we have 

(6-2) E (”) liBiBn-i{x) = ^{2x - l)S„_i(2a;) - ^S„(2 t) - Ji?„_i(cr). 
1^0 ^ ^ 

Proof. With a = b = e, the following limits are obvious: 

{e)i{e)n-i 1 2e(£)/ 1 1 

s-^o^ (2£)„ 2" E^oo (2e);_|_i 2* E->oo (2e -I- l)(2e) 4 

Hence we have 

(®-3) ^E (^i)Bi{x)Bn-i{x) =^2 (fi)fBiBr,-i{x) + ^B„_i{x). 

The sum on the left has a well-known evaluation (see (50.11.2) in [H]): 

E Bi{x)Bn-i{x) = n{2x- l)Bn-i{2x) - (n - l)H„(2a;). 

This, with (16.31) . immediately gives (16.21) . □ 

We note that (16.31) can also be obtained as a special case of identity (6.1) in [3]. 
For the next statement, again using a = b = e, we need the second-order harmonic 
numbers, defined by Fb^^^ := 0 and 

1=1 

We can now prove the following result. 
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Corollary 7. For integers n > 1 we have 


(6.4) = 


z=i 


Z=1 


HiF 


l\ B, 


I / I 


■Bn-l(x) 


+ nBn-ii^) + i Bnix). 


Proof. Setting a = 6 = e in (12.51) and dividing both sides by e, we have to expand 
the following terms. First, for 1 < Z < n — 1 we have 

(6.5) 

1 {o,)i{b)n-i (e + 1)... (e + Z — 1)(£ + 1)... (e + n — Z — 1) 


£ {a + b)r 


2{2s + l)...{2s + n-l) 


{I - l)\{n - I - ly. 
2(n- 1)! 

{i-iy.in-i-iy 




J 




i=i 


i=i 


2e 


2(n- 1)! 

Next, for Z > 1 we get 

1 a{b)i + b[a)i (e + 1)... (s + Z — 1) 1 


(l + {Hi-i + FLn-i-i — 2F[n-i) e + 0(£^)) . 


( 6 . 6 ) 


e (a + b)i^i (2e + 1)... (2£ + Z) Z 

z-i 


=i('+T) n 


-1 i-i 


1+7 


. , 1 + 2? 
1=1 1 




z 


2 £ 


1=1 


Z 


— - (l-(i7; + -)e + O(e^) ) , 


where we have used the fact that i7z_i + 2/1 = Hi + 1/Z. Next, we have 


(6.7) 


1 


ab 


1 


1 


= (l_2£ + 0(e2)). 


e (a + & + l)(a + 6) 2(1 + 2e) 2 

Finally, we collect on the right the terms left out in the two sums, namely 


1-2 


{a)r 


(a + b)r 


Bn{x). 


We expand 


, (a)n 2.(^s')r. 

(a + b)n (2£)r; 


n—1 


n 1+ 


i=i 


n—1 

n 

1=1 

e 

j 


£ + j 

2e + j 


1 + — 

J 


=n(i-+%-+-^ 

7=1 V 


— 1 - + ‘^Hn-l + X] “P 


l<j<k<n 
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Now the double sum in the last term can clearly be written as 
and thus 


( 6 . 8 ) 


1 - 2 


(a)^ 


(a + h)r, 


= . - 


2^n-l + 2^n-l ) £ + 0{e^). 


If we substitute (lOD-dOD into (I2.5p and let £ —0, we recover the polynomial 
analogue of Miki’s identity. Finally, if we equate the coefficients of e, we immediately 
get (16.41) after multiplying both sides by —1 and exploiting symmetry in (16.51) . □ 


6.2. Consequences of Theorem 2. We restrict our attention to the case fc = 3 
and oi = 02 = 03 = e. Furthermore, to avoid double indices, we set i = li, j = I 2 , 
I = I 3 on the left, and i = Iq, j = h, I = I 2 on the right of 611). Then we get, after 
dividing by n!, 

. . {e)i{e)j{e)i B,{x) Bj{x) Bi{x) _ £(£)j(£); B,{x) Bj Bi 

^ ^ (3s). ^! .■! d (3 £),+i+i 7! ,! /! 

,0 V £^(£)j Bi{x) Bj £3 Bn- 2 {x) 

(3s)i+ 2 *! j! ^(3£)3 (n-2)!’ 

valid for all n > 2 and £ > 0. This will be the basis for the various results in this 
subsection, and also immediately gives (63- 

For a first easy consequence we let £ — >■ 00 on both sides of (16.91) . Then the limit 
of the four fractions involving £ are easily seen to be 3“", 3“-’“^, and 3“^, 

respectively. Thus, after multiplying both sides by 3"n!, we immediately get the 
following result. 

Corollary 8. For integers n > 2 we have 
(6.10) ^ 

i+j+l=n \*jj! / 

.-1 / _ 1 \ 

(” . j3*Bi(x).B„_i_i + n(n - l)3”“^B„_2(a;). 

i=o k * / 

For the next consequence of (16.91) we set x = 0, for greater simplicity of the 
statement. The proof is tedious, and we leave the details to the interested reader. 


Corollary 9. For integers n > 2 we have 

1 
3 


1 Bi Bj Bi _ 

T 2-^ j j I 2-^ 


i-l) i j I ^ 


n-l\Bi Bn-i-i 


Y, +a - 2 x: () {^m +« 


n —1 


Z=1 
n — 1 




n — 1 


/ - 1 


/ + 1 y I n — I — 1 

Bl Bn-l 


I n — I 


-Bn_2 


1 


(n — l)n 


- 3 ) - 2 ( 


This can be seen as a third-order analogue of Miki’s identity. Note the difference 
in complexity between this result and the third-order analogue of Matiyasevich’s 
identity given in (14.21) . See also [3J (6.5)] for a third-order “Miki analogue” for 
Bernoulli polynomials. To prove Corollary 9, one can use a similar method as in 
the proof of Corollary 7, and proceed as follows: 
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- Collect the “edge” and “corner” terms in the sums of (I6.9|l . 

- Divide both sides of (j6.9l) by e^. 

- Expand the various fractions involving e in a similar way as in (j6.5l) and (16.61) . 

- Equate the constant terms (i.e., let £ —>• 0) to obtain a first identity. 

- Equate the coefficients of e to obtain a second identity. 

Interestingly, in this case the first identity turns out to be equivalent to (11.311 . Miki’s 
original identity. The second one is Corollary 9. 

6.3. Consequences of Theorems 3 and 4. Given the similarities between The¬ 
orems 1 and 3, it is clear that Euler analogues of Corollaries 2-7 could easily be 
derived; recall that an analogue of Corollary 1 is already stated following Theo¬ 
rem 3. Here we restrict ourselves to only a few more consequences; we also skip the 
proofs which are again similar to the proof of Corollary 7. 

Corollary 10. For integers n >2 we have 



, OCJ En-i{x) ^ ^Hn-i En{0) 
+ ^tln-2 -;-r 4--- 


n — 1 n — 1 n 


and for n > 1, 




r{2) \ Ei{0) I 3^(2)'\ ^»-l-l(0) 

'“V n + l-l I ^ n \ ^ "/n-tl 


Finally, to obtain third-order analogues of Miki’s identity, we start with (j5.26l) 
for fc = 3 and follow the outline described after Corollary 9. This leads to the 
following identities. 

Corollary 11. For integers n>2 we have 






and 


1 

3 
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Using the special values (j5.6ll . numerous identities involving Genocchi and/or 
Euler numbers could also be obtained. 


7. Final remarks 


1. By choosing different values of the parameters in Theorems 1-4, many more 
identities for Bernoulli and Euler numbers and polynomials could be obtained, some 
relatively simple, and others of increasing complexity. We have shown in this paper 
that the original identities of Miki and Matiyasevich and their various extensions 
are special cases of very general class of identities. 

2. This is not the first common extension of the identities of Miki and Matiya¬ 
sevich. In fact, using generating functions, Dunne and Schubert [5] recently proved 
the following result. 


Theorem 5 (Dunne and Schubert). For any integer n >2 and real p > 0 we have 


17 D ^!\ ^21 .B 2 n- 2 / _ r,p T{2n + 2p) ^ T(j) + l)T(p+l) 

(7.1) }^(20p(2n - 21), — ^^^ - 2S2„ ^ 

1=1 '■ ' 1=1 


2 r(p + i) 


(2n) 


1=1 ^ 2 


(2n)! ^ r(2p + ; + l) 

2n\ r(p -I- 2/)r(2p -I- 2n) 


T{2p + 2l + l) 


-B2lB2n-2 


The case p = 0 gives Miki’s identity, while for p = 1 we get 


n 

Y.B2lB2n -21 

1=1 



+ 2ni?2nj 


which is equivalent to Matiyasevich’s identity. 

The identity dzU) actually follows from Theorem 1 if we take a = b = p and 
a; = 0, then replace n by 2n and extract the end terms in the sums. After some 
manipulations we then get 


(7.2) ^(2Z)p(2n-20; 




B2I B2n -21 
’21 2n- 21 




2 r(p + i) 

(2n)! 


r^(p) 

(2n)! 

^ 2 n\ r(p 


i: 

1=1 ^ 2 


((2p)2n - 2 (p)2„) 

2l)T{2p + 2n) 
V{2p + 2l + l) 


B2lB2n-2 


Comparing (EH) with (j7.2ll shows that after some simplification we have 


r(p + o 

^ r{ 2 p + i + i) 


r(p) r(p-i-2n) 
Ti2p+1) ~ pr{2p + 2n) ■ 


This last identity can be proved independently, for instance by manipulating the 
integral representation of Euler’s beta function. 
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